Study of transmission and reflection from a disordered lasing medium 
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A numerical study of the statistics of transmission (t) and reflection (r) of quasi-particles from 
a one- dimensional disordered lasing or amplifying medium is presented. The amplification is intro- 
duced via a uniform imaginary part in the site energies in the disordered segment of the single-band 
tight binding model. It is shown that t is a non-self-averaging quantity. The cross-over length scale 
above which the amplification suppresses the transmittance is studied as a function of amplification 
strength. A new cross-over length scale is introduced in the regime of strong disorder and weak am- 
plification. The stationary distribution of the backscattered reflection coefficient is shown to differ 
qualitatively from the earlier analytical results obtained within the random phase approximation. 

PACS Numbers: 42.25.Bs, 71.55.Jv, 72.15.Rn, 05.40.+j 
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In recent years the study of wave propagation in an 
active random medium |I]-|lO|] , in the presence of ab- 
sorption or amplification, is being persued actively. The 
light propagation in an amplifying (lasing) medium has 
its implications for stimulated emission from random me- 
dia. In the Schrodinger equation, to describe the ab- 
sorption or amplification, one introduces the imaginary 
potentials. In that case the Hamiltonian becomes non- 
Hermitian and thus the particle number is not conserved. 
It should be noted that in quenched random systems with 
imaginary potentials the temporal coherence of the wave 
is preserved in spite of amplification or absorption which 
causes a particle non-conserving scattering process. 

Some new results have been obtained in this area. In a 
scattering problem, the particle experiences a mismatch 
from the real valued potential to the imaginary valued 
potential at the interface between the free region and 
the absorbing (or amplifying) medium, and hence it tries 
to avoid this region by enhanced back reflection. Thus a 
dual role is played by imaginary potentials as an absorber 
(or amplifier) and as a reflector. When the strength of 
the imaginary potential is increased beyond certain limit, 
both absorber and amplifying scatterer act as reflectors. 
Thus the reflection coefficient exhibits non-monotonic be- 
havior as a function of the absorption (amplification) 
strength. Using the duality relations it has been shown 
that amplification suppresses the transmittance in the 
large length (L) limit just as much as absorption does 
Pj . There exists a crossover length scale L c below which 
the amplification enhances the transmission and above 
which the amplification reduces the transmission which, 
in fact, vanishes exponentially in the L — > oo limit. In 
contrast, reflectance saturates to a finite value. More- 
over, absorption and amplification of same strength (i.e., 
differing only in the sign of the imaginary part) will in- 
duce same localization length |j] . 

In an amplifying medium even though the transmit- 
tance (i) decreases exponentially with the length L in the 
large L limit, the average (t) is shown Q| to be infinite due 
to the less probable resonant realizations corresponding 
to the non Gaussian tail of the distribution of In t. This 
result is based on the analysis using random phase ap- 



proximation (RPA). Using duality argument Paasschens 
et. al. show that non-Gaussian tails in the distribution 
of Int contain negligible weight |0j . Thus one might ex- 
pect finite value for (t) in the asymptotic limit. It should 
be noted that even in the ordered periodic system all 
the states are resonant states and still the transmittance 
decreases exponentially for all the states in the large L 
limit. The above simple case may indicate that in the 
asymptotic limit (i) is indeed finite. One of our objec- 
tives in this paper is to study the behavior of the trans- 
mission probability as a function of length in the presence 
of coherent amplification. We show that the transmis- 
sion coefficient is a non-self-averaging quantity. In the 
large length limit we do not find any resonant realization, 
which can give an enhanced transmission. We also study 
the behavior the of cross-over length L c as a function of 
amplification strength. We have analysed the behavior 
of logarithm of the transmittance which will have a max- 
imum value (lnt) max at L Cl and its dependence on the 
amplification strength. For a given strength of ampli- 
fication there exists a critical strength of disorder below 
which the average transmittance is always less than unity 
at all length scales and decreases monotonically. In this 
regime L c and (lnt} max loose their physical significance. 
In this regime we show that there exists a new cross- 
over length scale £ c which diverges as the amplification 
strength is reduced to zero for a given strength of the 
disorder. 

In the work by Pradhan and Kumar ||, the analyti- 
cal expression for the stationary distribution P s (r) of a 
coherently backscattered reflection coefficient (r) is ob- 
tained in the presence of amplification in a RPA using 
the method of invariant imbedding. The expression for 
P s (r) is given by 
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where D is proportional to rj/W, ij and W being the 
strength of amplifying potential and disorder respec- 
tively. One can readily notice from Eqn. ([!]) that P s (r) 
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does not tend to S(r — 1) in the large 77 limit. In this 
limit, as mentioned earlier, an amplifying scatterer acts 
as a reflector. The validity of above expression is limited 
to small disorder and amplification strength. 

We consider a quasi-particle moving on a lattice. The 
appropriate Hamiltonian in a tight-binding one-band 
model can be written as 

H = Y,e' n \n){n\ + V{\n){n + 1\ + \n)(n-l\). (2) 

V is the off-diagonal matrix element connecting nearest 
neighbors separated by a lattice spacing a (taken to be 
unity throughout) and \n) is the non-degenerate Wannier 
orbital associated with site n, where e' n = e„ — in is the 
site energy. The real part of the site energy e„ being ran- 
dom represents static disorder and e n at different sites are 
assumed to be uncorrelated random variables distributed 
uniformly (P(e„) = 1/W) over the range — W/2 to W/2. 
We have taken imaginary part of the site energy 77 to be 
spatially uniform positive variable signifying amplifica- 
tion. Since all the relevant energies can be scaled by V, 
we can set V to unity. The lasing medium consisting of 
N sites (n = 1 to N) is embedded in a perfect infinite 
lattice with all site energies taken to be zero. To cal- 
culate the transmission and reflection coefficients we use 
the well known transfer-matrix method , and the details 
are described in Ref. H . 

In our studies we have set the energy of the incident 
particle at E — 0, i.e., at a midband energy. Any other 
value for the incident energy does not affect the physics 
of the problem. In calculating average values in all cases 
we have taken 10,000 realizations of random site energies 
(e n ). The strength of the disorder and the amplifica- 
tion are scaled with respect to V, i.e., W (= W/V) and 
n (= rj/V). The length L = L/a. 

Depending on the parameters 77, W and L the trans- 
mission coefficient can be very large (of the order of 10 12 
or more). Hence, we first consider behavior of (Int) in- 
stead of (t). The angular brackets denote the ensem- 
ble average. The localization length £ can be computed 
from the behaviour of transmittancc in the large length 
limit. We denote the localization length i a for an or- 
dered medium (W=0)in the presence of uniform amplifi- 
cation. The localization length jllj for a disordered pas- 
sive medium (77 = 0) is given by elastic back scattering 
length I = 48V 2 /W 2 at the center of the band (E = 0). 
We have verified that the localization length in the pres- 
ence of both disorder and amplification, |[5j] ^ is related to 
I and l a (for n/V < 1 and W/V < 1) as £ = ll a /(l + l a ). 

In Fig. [l] we have plotted (lnt) max against 77 for a fixed 
value of W = 1.0 and the inset shows variation of L c with 
77 for W — 1.0. Initially (Int) increases with 77 and 
after exhibiting a maxima it decays to zero for large 77. 
This arises from the fact that the lasing medium acts as a 
reflector for large 77 as discussed in the introduction. Near 
the maximum, in a finite regime of 77, (lnt) max exhibits 
several oscillations. In this region sample to sample fluc- 
tuations of Int are very large. Thus average over 10, 000 



realizations may not represent the true ensemble aver- 
aged quantity. From the curve fitting of our numerical 
data for L c , we find that L c does not follow a power law, 
(1/^rj), in the full parameter regime ||. 

To study the nature of fluctuations in the transmis- 
sion coefficient, in Fig. || we have plotted, on log- 
scale, (t) , root-mean-squared variance t v — \J (t 2 ) — (t) 2 
and root-mean-squared relative variance (or fluctuation) 

t rv = J (t 2 ) — (t) 2 1 (t) as a function of L for 77 = 0.1 and 
W = 1.0. For these parameters I « 48,l a = 10, £ ps 8 
and L c ps 30. We notice that both (t) and t v exhibit 
maxima and decrease as we increase the length further. 
Except in the small length limit, variance is larger than 
the mean value. The relative variance is larger than one 
for L > 10 and remains large even in the large length 
limit. The t rv fluctuates between values 50 to 300 in the 
large length [L > 10) regime, indicating clearly the non- 
self-averaging nature of the transmittance. This implies 
that the transmission over the ensemble of macroscopi- 
cally identical samples dominates the ensemble average. 
In such a situation one has to consider the full probability 
distribution P(t) of t to describe the system behavior. 

We would now like to understand whether there ex- 
ist any resonant realizations in the large length limit for 
which the transmittance is very large. This study calls 
for sample to sample fluctuations. It is well known from 
the studies in passive random media that the ensemble 
fluctuation and the fluctuations for a given sample as a 
function of chemical potential or ener gy a re expected to 
be related by some sort of ergodicity |12| , i.e., the mea- 
sured fluctuations as a function of the control parameter 
are identical to the fluctuations observable by changing 
the impurity configurations. In Fig. ||(a) we have plot- 
ted t versus incident energy E (within the band from 
—2 to +2) for a given realization of random potential 
with n = and L = 100. The Fig. |(b) shows the 
behavior of t versus E for the same realization in the 
presence of amplification 7; = 0.1 and L = 100. ^From 
the Fig. [3](a) we observe that at several values of en- 
ergy the transmittance exhibits the resonant behavior in 
that t = 1. ^From Fig. ||(b) we notice that in the pres- 
ence of amplification, transmittance at almost resonant 
realizations is negligibly small. Few peaks appear in the 
transmittance whose origin lies in the combined effect of 
disorder and amplification. However, we notice that the 
transmittance at these peaks is much smaller, where as 
one would have naively expected the transmittance to be 
much much larger than unity in the amplifying medium. 
We have studied several realizations and found that none 
of them shows any resonant behavior where one can ob- 
serve the large transmittance. The peak value of ob- 
served transmittance is of the order of unity or less. This 
study clearly indicates that (t) is indeed finite contrary 
to the earlier predictions based on RPA |t| . 

So far our study was restricted to the parameter space 
of W and 77 for which L c and hence (lnt) max exist. 
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In Fig. [| we have plotted (Int) against L for ordered 
lasing medium (W = 0,77 = 0.01), disordered passive 
medium (W = 1.0, 7/ = 0) and disordered active medium 
(W = 1.0,r? = 0.01). The present study is restricted to 
the parameter space of 77 and W such that 77 <C 1.0 and 
W > 1.0. We notice that for an ordered lasing medium, 
the transmittance is larger than one. We have taken our 
range of L upto 300. For a disordered active medium 
(W — 1.0,7/ = 0.01), we notice that the transmittance 
is always less than one and monotonically decreasing. 
Initially, upto certain length, the average transmittance 
is, however, larger than that in the disordered passive 
medium (W — 1.0, 77 = 0). This arises due to the combi- 
nation of lasing with disorder. In the asymptotic regime 
transmittance of a lasing random medium falls below 
that in the passive medium with same disorder strength. 
This follows from the enhanced localization effect due to 
the presence of both disorder and amplification together, 
i.e., £ < I. It is clear from the figure that (Int) does 
not exhibit any maxima and hence the question of L c or 
(lnt) max does not arise. We notice, however, from the 
figure that for random active medium initially (Int) de- 
creases with a well defined slope and in the large length 
limit (Int) decreases with a different slope (correspond- 
ing to localization length £). Thus we can define a length 
scale £ c (as indicated in the figure) at which there is a 
cross-over from the initial slope to the asymptotic slope. 
In the inset of Fig. || we have shown the dependence of 
£ c on 77. Numerical fit shows that £ c scales as I/1/77, as 
we expect £ c — > 00 with 77 — > 0. As one decreases n, the 
absolute value of initial slope increases and that of the 
asymptotic one decreases. Simultaneously, the cross-over 
length £ c increases. In the 77 — > limit both initial as well 
as asymptotic slopes become identical. 

In Fig. ^ we have plotted the stationary distribution 
P s (r) of reflection coefficient r for different values of 77 (as 
shown in the figure) and a fixed value of W = 5.0. For 
small values of 77 = 0.05 the stationary distribution P s (r) 
has a single peak around r = r max — 1. The peak (r max ) 
shifts to higher side as we increase 77 (Fig. ||(b)). The 
behavior of P s (r) for small 77 is in qualitative agreement 
with Eqn. (Q). As we increase 77 further P s (r) exhibit a 
double peaked structure (Fig. §|(c)). At first the second 
peak appears at higher value of r at the expense of the 
distribution at the tail. As we increase 77 the second peak 
becomes more prominent and shifts towards left, where 
as the height of the first peak decreases. The distribu- 
tion at the tail has a negligible weight (see Fig. ||(c)). At 
still higher values of 77, the second peak approaches r«l 
whereas the first peak disappears. The now-obtained sin- 
gle peak distribution P s {r) in the large 77 limit tends to 
5(r — 1). In this limit the amplifying medium acts as 
a reflector and the disorder plays a sub-dominant role. 
The occurrence of the double peak structure along with 
P s (t) — > 5(r — 1) in the large 77 limit cannot be explained 
even qualitatively from Eqn. (|l|). This is due to the 
failure of RPA in this regime |3| . 

Our study for (lnr) 8 , obtained from P s (r) indicate that 



it initially increases with 77 and exhibiting a maxima at 
Vmax decreases monotonically. The double peak in P s (r) 
appears for values close to i] max . For the value 77 > rj max 
the amplifying medium acts dominantly as a reflector. In 
this regime we have verified that the stationary distribu- 
tion of the phase 8 of the complex reflection amplitude 
is not uniform. Physics of the double peak and overall 
shape of P s (r) shows similarity with the stationary dis- 
tribution obtained in the case of absorption (for details 
we refer to Ref. ||. 

In conclusion our numerical study on the statistics of 
transmission coefficient in random lasing medium indi- 
cates that in the asymptotic regime the transmission co- 
efficient is a non-self-averaging quantity, however, with 
a well defined finite average value. In some parameter 
space transmittance initially increases with 77 and falls 
off exponentially to zero in the asymptotic regime. In 
this regime there is a well defined cross-over length L c at 
which the transmittance is maximum, and it decreases 
monotonically with 77. In the parameter range where 
ij < 1, in the presence of disorder the average transmit- 
tance decreases monotonically and has a magnitude less 
than unity. In this regime L c does not exist. However, 
one can still define a new length scale £ c which scales as 
l/yfrj. Our study on the stationary distribution of re- 
flection coefficient P s (r) indicates that earlier analytical 
studies fail, even qualitatively, to explain the observed 
behavior in the large 77 limit. Our study clearly brings 
out the dual role played by an amplifying medium, as an 
amplifier as well as a reflector. 
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FIG. 2. The plot of (t), root-mean-squared variance (t v ) 
and root- mean-squared relative variance (t rv ) as a function 
of length L for tj = 0.1 and W = 1.0. 



FIG. 4. Variation of (Int) with L. The new length scale £ c 
which arises for r\ <C 1.0 is shown by a vertical dotted line. 
The inset shows the variation of £ c with r\ for W = 1.0. The 
numerical fit shown by the thick line indicates that £c scales 
as in this regime. 
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FIG. 5. Stationary distribution of reflection coefficient 
Ps(r) for W — 5.0 and various values of rj. The numerical 
fit shown in Fig. ^(b) with a thick line has D = 1.235 
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